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摘 要 : 为 探究 青海 省 植被 固 碳 量 时 空 演变 及 其 驱动 因子 ,基于 2000 一 2020 年 净 初 级 生产 力 (NPP) 数 据 和 土壤 呼 
吸 模型 计算 得 到 植被 净 生 态 系统 生产 力 (NEP) ,采用 趋势 分 析 、 偏 相关 分 析 及 地 理 探 测 吉 等 方法 对 青海 省 植被 NEP 
时 空 演变 及 驱动 因子 进行 定量 分 析 。 结 果 表 明 :青海 省 植被 NEP 在 近 20 a 间 波动 上 升 , 年 均 增幅 为 1.54 g C*m7*a^5 
年 均 植 被 NEP 空 间 变化 差异 较 大 , 由 东南 向 西北 递减 ,71.08% 的 区 域 保 持 不 变 或 增加 ;对 植被 NEP 解 释 力 最 强 的 是 
NDVI, 降 水 气温、 人 口 密度 等 气候 与 人 为 因子 对 NEP 的 空间 分 异 解释 力 较 强 ;由 于 双 因 子 交互 作用 会 增强 对 植被 


NEP 空间 分 异 的 解释 力 , 因 此 在 未 来 提升 青海 省 固 碳 能 力 时 ,要 注意 多 因子 协同 作用 。 
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陆地 生态 系统 是 全 球 气 候 变 化 研究 的 核心 内 
容 之 一 ,在 全 球 碳 平 衡 中 发 挥 着 重要 作用 '。 定 量 
估算 陆地 碳 汇 时 空 变 化 是 预测 气候 变化 和 执行 温 
室 气体 控制 协议 的 基础 。 通 过 比较 区 域 间 的 碳 乎 
衡 才 能 次 刻 理 解 区 域 间 循 环 过 程 。 通 过 人 工 生态 
工程 的 实施 来 提高 植物 固 碳 量 , 是 区 域内 增 汇 的 主 
要 手段 。 净 生态 系统 生产 力 (NEP) 是 单位 时 间 单 位 
面积 上 植物 的 碳 收文 状况 ,可 作为 定性 和 定量 化 分 
析 陆 地 植被 碳 循环 的 关键 参量 ,是 近 几 年 的 研究 
热点 ,如 基于 气候 变量 的 统计 模型 ,基于 光 能 利用 
率 相关 的 参数 模型 和 生态 系统 机 制 的 过 程 模型 等 ， 
尤其 是 3S 技 术 推 动 遥 感 反 演 获取 相关 地 表 参 数 和 
植被 信息 发 展 ,使 得 3S 技 术 支 撑 下 快速 发 展 的 生态 
系统 机 制 模型 成 为 了 当前 研究 的 主流 之 一 ”。 

青海 省 位 于 "第 三 极 " 东 北部 ,被 誉 为 “三江 之 
源 ,中 华 水 塔 ” ,是 对 全 球 气候 变化 响应 最 敏感 的 地 
区 之 一 , 故 青 海 省 植被 固 碳 能 力 变化 对 青藏 高 原 乃 
至 全 球 控 制 温室 气体 排放 和 碳 循环 均 有 巨大 影 
响 ”“。 在 青海 省 城市 化 和 工业 化 进程 中 ,人 口 分 布 、 
地 表 覆 盖 也 随 之 发 生 改 变 , 这 必 将 在 一 定 程度 上 影 
响 陆 地 上 的 植被 碳 循环 。 虽 然 ,已 确定 陆地 生态 
系统 将 在 未 来 持续 长 时 间 发 挥 碳 汇 的 作用 ,但 植被 
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NEP 驱动 机 制 和 植被 固 碳 量 变化 仍 存在 较 大 争议 ， 
如 何在 复杂 的 变化 环境 中 阐明 气候 敏感 区 植被 NEP 
近 几 十 年 的 变化 特征 以 及 驱动 机 制 有 竺 探究” 3。 刘 
凤 等 “和 周 夏 飞 等 对 青海 高 原 植被 青藏 高 原 草 
地 的 NEP 进 行 气 候 驱 动力 分 析 , 指 出 气候 变化 对 该 
区 域 NEP 的 演变 产生 影响 ; 张 新 中 等 ”通过 分 析 石 
羊 河流 域 的 NEP 与 NDVI. 降 水 .气温 和 太阳 辐射 的 
关系 ,指出 流域 内 NEP 变 化 主要 受降 水 和 和 气温 影 
响 , 受 太阳 辐射 的 影响 不 显著 ; 刘 春 雨 "研究 了 西 
南 地 区 植被 NEP 的 时 空 演 变 , 并 结合 人 类 强度 指 
数 , 采 用 地 理 加 权 回 归 法 探析 不 同 视角 下 的 人 为 因 
素 对 植被 NEP 影 响 。 已 有 对 植被 NEP 的 驱动 因素 
人 研究 主要 是 基于 气候 因素 ,对 人 为 因素 的 驱动 研究 
还 限于 单 因 子 的 影响 。 为 探析 双 因 素 交 互 作用 对 
植被 NEP 的 影响 ,本 文 将 地 理 探测 器 用 于 NEP 驱 动 
因素 的 定量 分 析 ,该 方法 既 可 进行 单 因子 影响 分 
析 , 亦 可 探测 多 种 人 为 因子 和 气候 因子 交互 作用 对 
NEP 的 影响 ”5 。 

近年 来 为 应 对 西北 地 区 “气候 转型 ”和 人 为 
活动 的 影响 ,研究 环境 变化 下 的 植被 固 碳 能 力 成 为 
实现 “ 矶 中 和 ”的 迫切 需求 ""。 本 文 利用 经 验 公 式 
和 MODIS 数 据 估算 青海 省 NEP, 并 探析 青海 省 NEP 
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时 空 演化 规律 及 驱动 因子 ,以 期 科学 评估 青海 省 植 
被 碳 储 量变 化 和 对 青海 省 植被 增 汇 政策 提供 理论 
参考 。 


1 研究 区 概况 与 数据 


1.1 研究 区 概况 

青海 省 位 于 89935'-103*04' E, 31?9' -39?19' N 
之 间 , 属 高 原 温 带 半 干旱 地 区 ,多 风 ,日 照 时 间 长 ， 
降雨 少 ,年 降水 量 在 50~550 mm, 年 均 温 在 -5.1~ 
9.0 C ;地 势 总 体 西 北向 东南 逐渐 倾斜 ;区 域 降 水 稀 
少 ,主要 集中 于 春 夏季 ,径流 主要 来 源 于 冰雪 融 水 
和 西南 暖 湿 气 流 。 

青海 省 林地 资源 主要 分 布 于 江河 源 区 附近 ,全 
省 林地 面积 2.65x10' km? ,天 然 草地 面积 4.2x10 
km 。 草 地 类 型 主要 为 高 寒 草 铅 和 草原 (图 1) ,主要 
分 布 在 祁连山 东北 部 和 上 某 达 木 盆地 东部 地 区 。 由 
于 特殊 的 气候 .高 海拔 和 风化 侵蚀 地 貌 ,导致 青海 
省 植被 易 受 人 类 和 和 气候 变化 影响 ”。 
1.2 数据 来 源 与 处 理 

净 初 级 生产 力 (NPP) 数 据 来 源 于 MODIS (https:// 
ladsweb.modaps.eosdis.nasa.gov/) 所 提供 的 MOD17A3- 
HGF 产品 ,时 间 分 辨 紊 为 年 ,空间 分 辨 率 500 m, 时 
间 跨 度 2000 一 2020 年 。 归 一 化 植被 指数 (NDVI) 数 
据 来 源 于 MOD13Q1 产品 ,时 间 分 辨 率 16 d, 空 间 分 
辩 率 250 m, 时 间 跨 度 2001 一 2020 年 ,采用 最 大 值 合 
成 月 数据 。 土 地 利用 数据 来 源 于 MCD12Q1 产品 ， 


时 间 分 辩 率 为 年 ,空间 分 辩 率 500 m, 时 间 跨 度 
2001 一 2019 年 。 以 上 数据 均 利 用 MRT 及 ArcGJS 进 
行 数据 拼接 ,并 投影 裁剪 研究 区 范围 。 

青海 省 气象 数据 来 自 2000 一 2020 年 国家 气象 
科学 数据 中 心 (http://data.cma.cn/) 月 值 数 据 集 ,空间 
分 辨 率 为 0.5°x0.5° ,对 该 数据 用 ArcGIS 进行 数据 还 
原 并 采用 样 条 也 数 法 插值 。 

人 均 GDP 和 人 口 密度 数据 均 来 自 中 国 科学 院 
资源 环境 科学 与 数据 中 心 (https://www. resdc.cn), 分 
辩 率 为 1km。 以 上 所 有 数据 均 采 用 Albers 等 面积 
圆锥 投影 进行 统一 分 析 。 

用 地 理 探测 需 探 测 空间 分 异性 时 ,选取 40 km 
的 空间 网 格 来 分 析 变 化 环境 对 植被 净 生 态 系统 生产 
力 (NEP) 的 影响 ,降水 .气温 采用 自然 间断 法 划分 为 
6 类 ;人 均 GCDP、 人 口 密度 采用 分 位 数 划分 为 6 类 ”。 


2 模型 与 方法 


2.1 NEP 的 估算 

NEP 常 被 用 来 描述 区 域 碳 平衡 动态 变化 ,可 指 
示 单 位 时 间 单 位 面积 上 植物 的 碳 收 支 状况 ,不仅 可 
作为 衡量 植被 活动 的 关键 因子 ,还 能 用 来 表征 植被 
的 健康 程度 和 固 碳 能 力 。 在 不 考虑 其 他 自然 及 
人 为 因素 影响 下 ,可 用 NEP 值 表征 陆地 和 大 气 生态 
系统 的 碳 交换 率 。NEP>0, 说 明 植 被 固 碳 量 多 于 土 
EIFE UR ( Ry ) 排 放 的 碳 量 ,表现 为 植被 碳 汇 ”Yi。 计 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2019)3333 号 ,对 底 图 边界 无 修改 。 下 同 。 


图 1 研究 区 植被 概况 
Fig. 1 Vegetation map of study area 
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NEP(x,t) = NPP(x,t) - R, (xt) (D — 据 研究 需要 ,本 文 应 用 因子 探测 器 ,交互 探测 器 等 4 


式 中 : NEP(x,t), NPP(x,t) FIR, xut) 为 第 1 年 像 元 x 
上 的 植被 净 生 态 系统 生产 力 .植被 净 初 级 生产 力 和 
土壤 呼吸 ,其 单位 均 为 gC.m。 

本 文 土壤 呼吸 尼采 用 裴 志 永 等 2 建立 的 经 验 
公式 ,由 气温 和 降水 间接 估算 。 经 已 有 研究 实地 采 
样 验证 ,该 公式 可 用 于 监测 青海 地 区 ,西北 地 区 士 
壤 呼 吸 。 计 算 公 式 如 下 3. 
R,=0.22x[EXP(0.0913x T) + In(0.3145 x P + 1)] x 


(2) 
30x 46.596 
式 中 :7 指 年 均 温 (%C);P 指 年 降水 量 (mm); Ry 为 年 土 
壤 呼 吸 (g C+m?)。 


2.2 驱动 因素 分 析 

22. 趋势 分 析 3H] Mann-Kendall 37" J£ Theil- 
Sen 法 分 析 2000 一 2020 年 植被 NEP 的 年 际 变化 规律 
与 动态 变化 ” ,计算 公式 如 下 : 


3 
var(S) 
Z6= 0, S-0 (3) 
— 
var(S) 
n-l nm 
S= b3 * sgn(x, — x;) (4) 


izlk-i*l 


XUP x, , ox; 为 连续 的 NEP 数 据 序列 ;为 年 份 。 标 
准 化 检验 统计 量 Ze ,其 中 Z .>0 时 表示 上 升 趋势 ， 
反之 , 则 表示 下 降 趋势 。 由 于 Mann-Kendall 法 强调 
单调 趋势 ,结果 不 一 定 是 线性 ,为 体现 变化 幅度 引 
和 Theil-Sen 斜率 法 表征 变化 幅度 Q KRN: 
Q=Median|(x, -x)/(k - i). Vk» i (5) 
式 中 :1<i<kp<n。 当 O >0 时 ,反映 上 升 趋势 ;反之 , 则 
反映 下 降 趋势 。 
2.2.2 相关 性 分 析 
(D) 偏 相关 系数 。 运 用 偏 相关 分 析 可 定量 刻画 
植被 NEP 对 其 他 影响 要 素 的 响应 。 偏 相关 系数 计 
P,-P,P, 
Pp: 5 (6) 
| (1- P) (1 - P?) 
假设 有 7 3 个 变量 , P;, 是 自 变 量 上 与 自 变 
量 i 和 j 的 偏 相关 系数 ; P; , Pao Pi 为 两 因子 间 的 
相关 系数 ”。 
(2) 地 理 探测 器 。 地 理 探测 器 是 一 种 探测 空间 
分 异性 ,并 揭示 其 背后 驱动 力 的 统计 学 方法 中 。 根 


种 探测 器 探测 NEP 的 驱动 力 及 多 因子 交互 作用 ，。 
子 蕊 对 属性 了 的 空间 分 异 解释 程度 用 4 值 度量 ,4 
值 的 计算 公式 如 下 : 


SN. h SSW 

LÉ4. $=! 三 十 三 

q-1 N, Le (7) 
L 

SSW= Y No, SST-N,, (8) 
h=1 


式 中 :9 为 某 因 子 对 NEP 空 间 分 异 的 解释 力 ,g 的 值 
域 [0,1],g 值 越 大 表征 该 因子 的 空间 分 异性 愈 明 
显 。 在 极端 情况 下 ,9g=1 表明 该 因子 完全 控制 了 
NEP 的 空间 分 布 。8 取 值 范围 为 [1 ,了 |] ,表征 因子 
的 分 区 , o, FIL o^ 2 3029 88 h ERAS XC NEP 的 方 
差 ,和 ,和 NN 分 别 为 层 h ME K HY f 762 ; SSWA SST 
分 别 为 层 内 方差 之 和 和 全 区 总 方差 。 
交互 探测 器 可 用 来 评估 不 同 因子 间 交 互 作用 对 
NEP 分 异 解释 力 的 强 弱 。 评 佑 方法 是 先 分 别 计算 2 
影响 因子 Xl1 和 了 对 NEP 的 解释 力 q, 再 计算 其 交 
互 作用 时 的 g 值 ,最 后 进行 比较 。 判 断 依据 如 表 1。 


表 1 交互 作用 判断 依据 


Tab.1 Basis of interaction judgment 


判断 依据 交互 作用 
q(X1 OX2)«Min[q(X1),4(42)] 非 线性 减弱 
Min[q(XD),a(X2)]< q(X1 (1X2) « 单 因 子 非 线性 减弱 
Max[q(X1).4(X2)] 
q(X1N X2) > Max[q( X1), q(X2)] 双 因 子 增强 
q(X11X2) 2 q(X1) + q(X2) 独立 
q(X11X2) » q(X1) + q(32) 非 线性 增强 


3 结果 与 分 析 


3.4 NEP 精 度 评价 

由 于 NEP 的 实测 数据 难以 获取 ,本 文采 用 其 他 
学 者 在 青海 省 实测 NEP 结 果 .青藏 高 原 数据 中 心 三 
江 源 气 象 站 点 实测 NPP 结 果 和 站 点 气象 数据 模拟 
Rn 结果 计算 NEP””。 将 实测 数据 与 估算 结果 进行 
相关 性 分 析 、 显 著 性 检验 ,表明 两 者 的 相关 系数 为 
0.864( 图 2) ,实测 与 估算 数据 显著 相关 (R=0.7821， 
P<0.01) ,显示 结果 具有 可 靠 性 。 
3.2 青海 省 NEP 时 空 分 布 特征 
3.2.1 青海 省 NEP 空间 分 布 特 征 ” 近 20 a 青海 省 年 
3 NEP 的 空间 分 布 差异 显著 ,呈现 出 东南 高 西北 低 
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2 青海 省 NEP 模 拟 结果 验证 
Fig. 2 Verification of NEP simulation results in 


Qinghai Province 


的 特征 (图 3)。 青 海 省 年 均 植 被 NEP 为 134.05 g C- 
ma', 其 中 NEP>0 区 域 高 达 4.89x10;km( 占 比 
74.15% ) ,NEP<0 区 域 可 达 890 km。NEP 高 值 区 
(312~735.90 g Cm”a ) 面 积 总 占 比 5.07% ,主要 位 
于 祁连山 附近 青海 湖 湿地 附近 、 江 河源 区 和 黄河 
上 游 谷地 等 地 。 植 被 碳 源 区 主要 有 格尔木 市 . 治 
县 和 杂 多 县 等 地 。 植 被 在 空间 分 布 上 过 渡 明 显 的 
原因 ,主要 与 当地 的 降水 状况 以 及 植被 分 布 有 关 。 
NPP 多 年 均值 分 布 与 NEP 分 布 相似 。 应 呈 西 北向 
东南 递增 之 势 , 具 体 分 布 为 森林 > 草原 > 荒漠 ,空间 
差异 较 明 显 。 

3.2.2 青海 省 NEP 的 年 际 变化 特征 ” 近 20a 青 海 省 
NEP 年 际 变化 较 明 显 (图 4) ,2000 一 2020 年 青海 省 
生态 系统 的 固 碳 能 力 呈 波动 增加 趋势 。 年 均 植 被 
NEP 1E 114.55-155.49 g C- m?-a ,NEP 年 均 增 幅 为 
1.54 gC.m2.a:。 植 被 NEP 的 极 大 值 、. 极 小 值 分 别 
出 现在 2018 年 .2001 年。 按照 植被 NEP 近 20 a 的 变 


(a) (b) 


mz 


NEP 年 均值 /(g C-m?-a7) 
O 无 植被 区 


0 200km - 735.90 0  200km 
-15.02 


Ryl(g Cm2.a) 


21.32 
- 8.35 


化 态势 ,划分 3 个 时 段 讨论 其 年 际 变化 特征 :2000 一 
2008 年 间 植被 NEP 波 动 增加 ,2006 年 出 现 第 1 个 峰 
值 144.93 g C*m?-a;2008—2015 4N] , NEP HZ 
年 均值 上 下 波动 ;2015 一 2020 年 间 ,植被 NEP 值 逐 
年 上 升 ,该 阶段 也 为 增幅 最 大 的 时 段 。 
3.2.3 青海 省 NEP 的 动态 交 化 ”为 研究 青海 省 NEP 
的 空间 动态 变化 规律 ,采用 逐 像 元 Theil-Sen 斜率 人 
计 法 ,得 出 青海 省 21 a NEP 变 化 趋势 (图 5)。 植 被 
NEP 空间 分 布 增长 格局 与 其 空间 分 布 格局 相似 ,从 
变化 趋势 可 知 , 青 海 省 植被 固 碳 能 力 整 体 上 有 增长 
趋势 , 少 部 分 呈 负 增长 趋势 的 区 域 主要 分 布 在 海 西 
市 和 玉树 市 内 , 约 占 青海 省 面积 的 2.17% , 主要 为 植 
被 稀 琉 地 区 ;增长 幅度 大 的 高 值 区 ( 占 比 3.17% ) 主 
要 位 于 青海 湖 湿 地 ` 贵 南 县 和 共和 县 附近 等 地 ,这 
些 区 域 的 降水 较 多 ,植被 覆盖 丰富 。 将 Theil-Sen i 
势 分 析 和 M-K 显 著 性 检验 结果 霉 加 ,通过 P<0.05 显 
著 性 水 平 的 区 域 面积 为 3.55x10; km ,其 中 显著 增加 
区 域 占 全 省 面积 的 $2.95% , 主要 分 布 在 柴 达 木 贫 地 
荒漠 ,祁连山 森林 与 高 寒 草原 .江河 源 区 -甘南 高 寒 
草 旬 草 原 等 生态 区 ;显著 减少 的 区 域 仅 占 0.23% , 零 
星 分 布 于 天 上 峻 县 BZH . 曲 麻 莱 县 EPK ` 治 多 
县 等 地 。 
3.3 青海 省 NEP 驱动 机 制 分 析 
3.3.1 NEP 对 气候 因子 的 响应 ”气候 因子 通过 影响 
植被 的 生长 态势 . 净 初 级 生产 力 以 及 土壤 呼吸 , 进 

影响 植被 净 生 态 系统 生产 力 的 变化 ;尤其 是 气 
温 、 降 水 这 2 个 主要 气象 因素 ,在 一 定 程度 会 影响 植 
被 NEP 的 变化 。 

由 图 6 可 知 , 近 20 a 青海 省 植被 NEP 的 极 大 值 

出 现在 2010 年 ,降水 量 极 大 值 出 现在 2009 年 


(c) 


mz 


NPP 均 值 /(g C-m?a7) 
O 无 植被 区 


0  200km 75521 
0.87 


图 3 青海 省 NEP .Ry 和 NPP 空间 分 布 
Fig.3 Spatial distribution of NEP, Ru and NPP in Qinghai Province 
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VÉ TIMES — 关 小 于 气温 的 偏 相关 。 初 步 分 析 王 旱 和 半 干 旱 区 
155 上 一 NEP 线 性 拟 合 线 。 有 iege : Sop 
iso L = 95% 置 信 区 间 气温 对 植被 的 分 布 及 长 势 影响 更 大 。 

$ 145 上 2000—2020 年 ,青海 省 降水 介 于 325-460.85 

CENAS a mm ,年 平均 降水 量 为 395.26 mm。 降 水 与 植被 NEP 

E 130 b 的 偏 相 关系 数 介 于 -0.24~0.42 之 间 ( 图 7a)。 植 被 

d NEP 与 降水 呈 负 相关 的 区 域 占 全 省 面积 的 23.32%， 
usL 主要 分 布 在 玉树 EHE RR, HEAR LIA IR . 班 玛 、 


110 


图 4 青海 省 NEP 的 年 际 变化 


Fig.4 Interannual variation of NEP in Qinghai Province 


(463.54 mm) ,在 2010 年 左右 年 均 温 也 有 1 个 显著 的 
峰值 。 说 明 青 海 省 在 2010 年 前 后 水 热 组 合 相 对 较 
好 ,有 利于 提高 生态 系统 固 碳 能 力 。NEP 的 极 小 值 
出 现在 2001 年 ,该 年 降水 量 显 著 偏 低 ,气温 虽 相 比 
2000 年 略 有 升 高 ,但 长 历时 干旱 易 造 成 土壤 呼吸 增 
强 和 植被 退化 ,从 而 导致 固 碳 能 力 下 降 。 对 上 述 3 
个 因素 间 的 响应 程度 进行 分 析 ,2000 一 2014 年 间 植 
被 NEP 变 化 相对 降水 存在 一 定 的 沾 后 性 ,滞后 期 约 
为 1 ai;2014 年 后 两 者 变化 相对 一 致 。 近 20 a ,气温 
与 植被 NEP 的 变化 保持 一 致 。 上 述 分 析 表 明 ,青海 
省 植被 NEP 对 年 均 温 和 年 降水 的 响应 程度 较 高 。 
空间 上 分 析 青 海 省 植被 NEP 与 年 均 温 、 年 降水 
量 的 偏 相 关 关 系 , 并 对 其 进行 P<0.05 的 显著 性 检验 
(图 7、 图 8)。 区 域内 植被 NEP 与 降水 的 偏 相 关系 数 
均值 为 0.134, 与 气温 的 偏 相 关系 数 均 值 为 0.429, 这 
表明 气温 降水 与 NEP 均 呈正 相关 , 且 降 水 的 偏 相 


AZ 河南 县 等 地 。 植 被 NEP 与 年 降水 量 呈 正 相关 
的 区 域 占 青海 省 面积 的 65.75% , 主要 分 布 在 海南 
市 海北 市 海 西 , 海 东 等 地 ;其 中 与 降水 呈现 显著 
正 相 关 (P<0.05 ) 的 区 域 面积 达 7.81 x 10* km’, 4] i 
全 省 面积 的 11.7%( 图 8a)。 

2000—2020 年 ,青海 省 年 均 温 -2.35~-0.70 C, 
年 均 温 均值 -1.44 %。 植 被 NEP 与 气温 呈 负 相关 的 区 
域 面积 为 1.142 x 10'km2, 仅 占 总 面积 的 1.71% ,零星 
分 布 于 各 市 ;植被 NEP 与 年 均 温 呈正 相关 的 区 域 占 
青海 省 面积 的 71.5% ,这 些 区 域 位 于 山药 处 .江河 源 
区 以 及 河流 谷地 ,主要 为 森林 ,草原 。 其 中 ,通过 显著 
性 检验 (P<0.05) 的 植被 NEP 与 气温 正 相 关 的 区 域 面 
积 达 2.572 x 10 km , 占 全 省 面积 的 38.53%( 图 8b)。 
3.3.2 _ NEP 与 NDVI 的 相关 性 ”分析 表明 ,青海 省 年 
均 植被 NEP 与 NDVI 的 空间 分 布 具有 高 度 相似 性 ， 
在 青海 湖 湿地 、 贵 南 且 和 共和 县 等 地 是 植被 NEP 
(312-735.90 g C- m^- a) Ej NDVI(0.68-0.91) HJ rej 
值 分 布 区 。 在 空间 分 布 上 (图 9) , TE NEP 5j NDVI 
呈正 相关 的 区 域 占 全 省 面积 的 65.83% ,其 中 通过 显 
著 性 检验 (P<0.05) 的 区 域 占 比 达 31.43% ,这 些 区 域 
主要 分 布 在 海南 、 海 东 和 海 西 东南 部 。NEP 与 NDVI 


基于 Theil-Sen 的 NEP 变 化 趋势/(g Cm2.aD 
mE 一 19~0 (3 1.1~2.4 mm 4.8-14 
æ 0.01-1 Eu 2.5-4.7 


通过 显著 性 95% 的 M-K+Theil-Sen 趋 势 /(g Coma) 
mm —19~0 C 1.1~2.4 mm 4.8~14 
æ 0.01-1 Eu 2.5~4.7 


图 5 青海 省 Mann-Kendall fil Theil-Sen fii NEP 的 空间 分 布 
Fig. 5 Spatial distribution of vegetation NEP of Mann-Kendall and Theil-Sen in Qinghai Province 
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-- 降水 量 -v- NEP 一 - 气温 m df FH É H9 pda 省 面积 的 7.37% , 主要 好 
» Ro med uc cer m 0 或 占 全 省 is isa es ciiin 
px De y=1.5428x+117.72 y-0.0349x-1.8325 UBER, TERI RE LA ih Tf DUAE S 
2—0.34 R2—0.5574 R?=0.3088 " 0.05 显著 1 竹 检 验 的 区 域 仅 占 0.1296, 这 些 区 域 零星 Ed 
500 "1 
4 168 分 布 于 各 地 。 
is (56 个 333 NEP 的 人 为 驱动 力 土地 利用 类 型 转变 能 揭 
È 300 & 1 144 号 ”示人 为 因素 对 植被 NEP 的 影响 。 参 考 《 土 地 利用 现 
$ 20 R o 9 状 分 类 ?和 原 数据 ICGBP 标 准 ,将 MCD12Q1 重 新 划分 
100 Lig 2. ， 为 6 大 类 (耕地 、 林 地 草地、 水域, 城乡 居民 用 地 ,未 
1 96 利用 地 )。 为 研究 2000 一 2020 年 土地 利用 类 型 对 青 
4 84 


2000 2005 2010 2015 2020 
年 份 


图 6 2000 一 2020 年 青海 省 NEP 与 年 均 温 、 
年 降水 量 的 年 际 变化 
Fig.6 Interannual variations of NEP, annual temperature and 


annual precipitation in Qinghai Province from 2000 to 2020 


降水 偏 相 关系 数 
mm —0.816- —0.218 m 0.202-0.427 
0 100km = —0.217-0 mm 0.428-0.879 
器 0.001-0.000 O 无 植被 区 


海 省 NEP 的 影响 ,结合 数据 的 时 间 跨 度 (2005 ^F. 
2010 年 .2015 年 .2020 年 ) ,以 Sa 为 1 个 研究 期 , 探 
究 人 为 因素 对 植被 NEP 的 影响 。 分 析 表 明 ,2015 一 
2020 年 的 土地 利用 类 型 转换 最 为 明显 ,转换 面积 达 
研究 区 总 面积 的 2.7% ,该 阶段 主要 是 未 利用 地 和 草 
地 之 间 的 类 型 转换 ,主体 为 荒漠 向 草地 转换 ,其 面 


人 


气温 偏 相关 系数 

EH _0.659~0 EU 0.425~0.571 
0 100km = 0.001-0.265 mm 0.572~0.901 
—= c 0.266-0.424 O 无 植被 区 


图 7 NEP 与 降水 和 年 均 温 偏 相关 关系 


Fig. 7 Correlation of NEP with precipitation and annual mean temperature deviation 


降水 偏 相 关 显 著 性 


0 100km EE 显著 正 相 关 显著 负 相 关 
i CO 不 显著 正 相 关上 国 不 显著 负 相关 


气温 偏 相 关 显 著 性 
0 100km EE 显著 正 相 关 


一 


E 显著 负 相 关 
CU 不 显著 正 相关 mb 不 显著 负 相关 


图 8 NEP 与 降水 和 年 均 温 偏 相 关 显著 性 分 析 


Fig. 8 Correlation analysis of NEP with precipitation and annual mean temperature deviation 
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(a) 


人 


NDVI 的 相关 系数 
EE —0.869-0 Œ 0.327~0.588 
E 0.001-0.326 mE 0.589-0.988 


0 100km 
L1 


NDVI 通 过 95% 显 著 性 水 平 的 相关 系数 
0 100km = —0.869— — 0.447 0.001—0.726 
= E -0.446-0 EH 0.727~0.988 


图 9 植被 NEP 与 植被 NDVI 相 关 性 系数 分 布 


Fig.9 Distribution of correlation coefficient between vegetation NEP and NDVI 


积 占 本 阶段 转换 总 面积 的 49.48% , 主要 位 于 柴 达 木 
盆地 荒漠 区 .昆仑 山东 段 高 寒 荒漠 草原 JU RE IR 
寒 荒漠 草原 区 , 主要 为 荒漠 与 草原 过 渡 区 域 。 耕 地 
面积 在 2005 一 2015 年 星 下 降 趋 势 并 转换 为 草地 , E 
要 是 政府 在 该 时 期 内 采取 了 退耕 还 草 等 生态 保护 
措施 。2001 一 2020 年 荒漠 与 草地 间 的 转换 较为 频 
繁 ,2010 一 2015 年 荡 漠 化 面积 略 大 于 绿洲 化 面积 ,该 
阶段 存在 微弱 的 荒漠 化 趋势 ,植被 NEP 呈 下 降 趋 势 
(图 10)。 


3.3.4 因子 影响 力 探测 “为 探究 不 同 因子 对 青海 省 
植被 NEP 的 影响 ,选取 6 种 典型 因子 并 对 其 进行 地 
理 探 测 器 分 析 ,探测 其 空间 分 异 的 驱动 力 ( 其 中 
2020 年 人 均 CDP 数 据 缺 失 )。 单 因子 对 植被 NEP 的 
影响 分 析 结 果 表明 (图 11) ,在 该 研究 期 内 ,NDVI 是 
影响 植被 NEP 的 主导 因子 ,其 q 值 均 在 0.7 以 上 ,并 
且 年 际 波动 较 小 ,说 明 植 被 覆盖 度 对 植被 NEP 贡献 
较 大 。 和 气温 与 降水 的 解释 力 较 强 , 均 超 过 了 0.3 , 降 
水 在 2015 年 的 影响 力 超过 和 气温, 主要 是 由 于 年 降水 


2005 年 2010 年 


2015 年 


2020 年 


mm 草地 mm 林地 c3 耕地 mm 未 利用 地 国 水 域 国 城乡 居民 用 地 


图 10 青海 省 土地 利用 转移 桑 基 图 


Fig. 10 Sangji map of land use transfer in Qinghai Province 
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降水 气温 密度 利用 NDVI 0.832 


0.787 0.831” 0.751' 0.746” 0.735 


0.693 


Hé 035* 0.137 


0.554 
0.414 
0.275 


2020 年 0.136 


注 :* 表 示 在 0.05 水 平 上 ( 双 侧 ) 显 著 相 关 。 


11 代表 年 份 植被 NEP 单 因子 影响 力 和 交互 作用 gq 值 


Fig. 11 Single factor influence and interaction q value of vegetation NEP in representative years 


偏 少 ,限制 了 植被 生长 ,从 而 抑制 了 植被 的 固 碳 能 
力 。 土地 利用 和 人 均 GCDP 影 响 力 相 对 稳定 且 解 释 
力 较 低 , 人 口 密度 解释 力 较 高 ,在 2010 年 的 影响 力 
高 于 气温 ;降水 在 2020 年 解释 力 略 高 于 土地 利用 。 

地 理 探测 器 交互 作用 结果 表明 (图 11): 区 域内 
两 因子 交互 作用 的 g 值 显著 大 于 单 因 子 的 g 值 ,表明 
交互 作用 呈 非 线性 增强 ,两 因子 交互 作用 后 对 植被 
的 NEP 解 释 力 会 显著 增强 。 对 4 个 代表 年 份 的 主导 
交互 因子 进行 排列 ,最 大 主导 作用 的 交互 因子 是 该 
年 最 强 主导 因子 与 其 他 因子 的 组 合 。2000 一 2020 
年 最 强 主导 交互 因子 均 为 NDVI 与 人 均 GDP 气温 
与 NDVI, 解 释 力 均 达 到 了 0.8 以 上 。 虽 然 CDP 的 单 
因子 影响 力 小 ,但 当 人 均 GDP 与 NDVI 交 互 作用 时 ， 
在 特定 年 份 会 超越 其 他 因素 。 


4 讨论 


MOD17A3HGF 模型 已 被 用 于 估算 青海 省 NPP, 
通过 作物 产量 验证 ,证 实 估算 结果 可 靠 " 。 估 算 
Rs 的 检验 模型 在 青海 高 原 和 青藏 高 原 通过 了 实地 
KUE 1。 

本 文 估算 青海 省 近 20 a NEP 年 均值 为 134.69 g 
Cm2.a ,与 周 夏 飞 等 (120.8 g Cem- a") 、 刘 风 等 ” 


(128.40  C- m^: a 0) 的 估算 结果 较为 接近 。 虽 然 本 
文采 用 的 NEP 估 算 方 法 不 同 且 更 新 了 时 间 序 列 ,但 
估算 结果 与 前 人 的 较为 接近 ,说明 佑 算 方法 可 行 。 
该 方法 存在 的 不 确定 性 ,主要 源 自 土壤 呼吸 的 不 确 
定性 。 将 应 定义 为 与 气象 因子 相关 的 函数 ,容易 忽 
视 呼 吸 底 物 的 分 解 对 异 养 呼吸 的 贡献 ,容易 造成 高 
植被 覆盖 度 区 域 的 碳 汇 高 估 , 低 植被 覆盖 度 区 域 对 
碳 汇 的 低 佑 ,因此 需 慎 重 采用 ”。 此 外 ,Rs 究竟 是 
否 适用 于 青海 省 全 省 植被 NEP 的 估算 需要 大 量 的 
地 表 观 测 实 测 值 作 参考 ,尽管 以 多 年 平均 值 做 比 
较 , 但 时 间 跨 度 的 区 别 难免 产生 一 定 误差 ”，。 本 文 
NEP 与 气温 的 关系 ,与 已 有 人 研究 存在 一 定 差异 ”, 主 
要 由 于 两 者 对 NPP 的 估算 方式 及 采用 的 气象 原 数 
据 插 值 方 法 存在 差异 所 致 ,在 中 亚 地 区 和 石 羊 河流 
域 的 研究 结论 与 本 文 相似 。 

诸 因 子 对 植被 NEP 的 影响 分 析 表 明 ,影响 较 大 
的 是 以 气温 、 降 水 为 代表 的 气候 因子 和 以 土地 利用 
为 代表 的 人 为 因子 。 分 析 结 果 显 示 ,气温 对 NEP 的 
影响 程度 要 大 于 降水 ,原因 是 适度 的 升温 对 NPP 的 
促进 作用 要 大 于 应 的 促进 作用 王 ,而 且 当 海拔 超 
过 2000 m 时 ,水 分 对 于 植物 NPP 的 限制 作用 越 来 
越 弱 交 。 与 气温 呈 偏 负 相 关 的 主要 是 荒漠 区 域 ( 如 


cr 
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柴 达 木 贫 地 荒漠 区 ) ,该 区 域 地 表 缺 水 .蒸发 旺盛， 
气温 降低 会 导致 区 域内 莹 发 减弱 ,有 利于 缺 水 状况 
缓解 并 促进 植被 NEP 的 积累 ,进而 提升 该 地 区 的 植 
被 固 碳 能 力 。 受 三 江 源 头 太阳 辐射 强烈 .冰川 资源 
丰富 ,海拔 高 等 因素 影响 ,在 青海 南部 NEP 与 降水 
呈 负 相关 .与 气温 呈正 相关 。NDVI 与 植被 NEP 在 
空间 分 布 上 具有 极 大 相似 性 ,植被 覆盖 率 高 的 区 域 
光合 效率 高 ,植物 群落 固定 在 地 面 的 有 机 碳 含量 及 
植被 净 碳 交换 量 同 样 也 高 。 

运用 地 理 探测 器 探测 青海 省 NEP 的 人 为 及 气 
候 驱 动力 ,整体 上 气候 对 植被 NEP 地 理 分 异性 的 解 
释 力 要 超过 人 为 因素 ,但 在 某 些 特殊 年 份 ,情况 过 
异 。 在 气候 趋 于 “ 暖 湿 ”背景 下 "9 ,由 于 气温 逐年 上 
升 ,陆地 生态 系统 的 植被 储 碳 能 力 对 气候 的 响应 将 
呈现 新 的 格局 ,人 为 因素 对 NEP 空 间 分 异 解释 力 较 
3g ,但 与 其 他 因子 交互 时 均 超 过 0.3。 双 因子 的 协同 
交互 作用 远 超 单 因 子 的 作用 ,未 来 在 提升 生态 系统 
固 碳 能 力 时 , 需 在 关注 当前 气候 变化 的 同时 ,采取 
适当 的 人 为 干预 。 降 水 对 NEP 的 解释 力 在 2020 4E 
最 低 的 原因 尚 待 进一步 探究 。 

由 于 生长 季 内 农作物 固定 的 碳 将 迅速 进入 大 
气 碳 循环 ,鉴于 研究 区 内 农田 比例 较 大 , 若 能 剔除 
作物 NEP 的 影响 将 对 区 域 增 汇 更 有 指导 意义 7。 
青海 省 作为 五 大 牧区 之 一 ,在 后 续 研 究 中 还 应 将 牲 
畜 饲 养 量 纳 人 人 为 因素 ,分 析 其 对 青海 植被 NEP 的 


影响 。 


5 结论 


“l 


本 文 详细 分 析 了 青海 省 近 20 a NEP 的 时 空 变 
化 规律 及 驱动 因素 ,并 利用 地 理 探测 器 定量 分 析 不 
同 因 子 对 NEP 的 影响 ,具体 结论 如 下 : 

(1) 时 间 尺 度 上 2000—2020 年 青海 省 NEP 呈 
增加 趋势 ,增幅 为 1.07 g Cema, SEFE NEP 为 
134.69 g Cm -a'o Z EE , NEP Œ A ZR I] P8 335 7 
的 特征 ,NEP 分 布 与 植被 覆盖 基本 一 致 。 

(2) 在 空间 与 年 际 尺度 上 ,植被 NEP 与 气温 、 降 
水 均 呈 正 相关 ,气温 与 NEP 的 相关 性 大 于 降水 。 

(3) 地 理 探测 器 表明 , 近 20 a 植被 覆盖 度 对 
NEP 起 主导 作用 ,气温 的 影响 作用 相对 较 强 ,在 干 
旱 年 份 降水 的 影响 力 要 超过 除 NDVI 以 外 的 其 他 因 
子 , 双 因子 的 交互 作用 超过 单 因子 影响 力 。 青 海 省 
植被 NEP 空间 分 布 是 气候 和 人 为 因素 共同 作用 的 
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结果 , 且 不 同 影响 因素 对 其 作用 的 强 弱 不 同 ,在 特 
定年 份 人 类 活动 与 其 他 因素 交互 作用 是 植被 NEP 
变化 的 主导 因素 。 
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Abstract: Based on NPP data from 2000 to 2020 and the vegetation net ecosystem productivity (NEP) calculated 
by the soil respiration model, the spatial-temporal adaptation and driving factors of vegetation NEP in Qinghai 
Province were quantitatively analyzed via trend analysis, 6 partial correlation analysis, and geographical detector 
to explore the spatiotemporal adaptation of vegetation carbon sequestration and its driving factors. The results 
showed that the vegetation NEP fluctuated over the past 20 years, with an average annual increase of 1.54 g C +m” 
a '. The spatial variation of annual vegetation NEP varied greatly, decreasing from southeast to northwest, and 
711.0896 of the area either remained unchanged or increased. Normalized difference vegetation index (NDVI) has 
the strongest explanatory power for vegetation NEP, and climate and human factors, such as precipitation, 
temperature, and population density are stronger factors for the spatial differentiation of NEP. Because the two- 
factor interaction will increase the strength of the argument for vegetation NEP spatial differentiation, it is necessary 
to pay attention to multi-factor cooperation in the future to enhance the sequestration capacity of carbon in Qinghai 
Province. 
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